Detecting unambiguously non-Abelian geometric phases with trapped ions 
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We propose for the first time an experimentally feasible scheme to disclose the noncommutative 
effects induced by a light-induced non-Abelian gauge structure with trapped ions. Under an ap- 
propriate configuration, a true non-Abelian gauge potential naturally arises in connection with the 
geometric phase associated with two degenerated dark states in a four-state atomic system inter- 
acting with three pulsed laser fields. We show that the population in atomic state at the end of a 
composed path formed by two closed loops C\ and C2 in the parameter space can be significantly 
different from the composed counter-ordered path. This population difference is directly induced by 
the noncommutative feature of non-Abelian geometric phases and can be detected unambiguously 
with current technology. 
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Gauge structures have been used to describe almost 
all the fundamental interactions in nature and have been 
shown to play a key role in modern physics. Gauge 
fields are usually classified as Abelian and non-Abelian 
ones according to the commutative rule of the associated 
group elements. Remarkably, they also arise naturally 
in the adiabatic evolution of simple quantum systems, 
whose initial formulation has no apparent relationship 
to gauge fields 1, 2|. Berry demonstrated that the adi- 
abatic cyclic evolution of a non-degenerated state leads 
to a scalar geometric phase factors corresponding to a 
U(l) Abelian gauge field Wilczek and Zee further 
generalized this idea to the case of degenerate states and 
showed that a non-Abelian structure emerges if a set of 
quantum states remains degenerate as the Hamiltonian 
varies. When a system with A^-fold degenerate levels un- 
dergos cyclic process, the mapping corresponds to a geo- 
metric phase factor with SU (N) x U(l) matrix formJJQ. 
So far, gauge structures in such simple quantum sys- 
tems have found man y a pplications in diverse fields in 

physicsi ansa sm eh, ee ee ei m . 

Following the progress achieved in theory and appli- 
cation, direct experimental observations of gauge struc- 
tures in quantum systems attracts increasing interests. 
Although Abelian geometric phases have been tested ex- 
perimentally by various means 

aa, 

the noncommuta- 
tive effect, a key property of non-Abelian structures has 
not been experimentally demonstrated. In an interest- 
ing experiment fl6j. Tycko demonstrated the effect of ge- 
ometric phase on the magnetic resonance spectrum of 
two-fold degenerated statesJUj]. However, as pointed out 
in Ref. 0], because of the experimental limitation in 
nuclear magnetic resonance (NMR), only the (commu- 
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tative) Abelian part of the non-Abelian gauge structure 
has been experimentally observed and the direct obser- 
vation of noncommutativity of non-Abelian gauge struc- 
tures in NMR system is really difficult. Recent exper- 
iments in trapped-ion systems provide high-precise ma- 
nipulations and measurements on internal states of atoms 
HE EE EE HE HE HE E3|. Comparing with NMR sys- 
tems, true non-Abelian structures can be easier to be re- 
alized by modulating controlling parameters in principle. 
Thus it is significant and important to explore whether 
a trapped-ion system is promising for detecting the fun- 
damental non-Abelian characteristic. 

In this paper, we analyze in detail why non-Abelian 
structures have not been observed in NMR experiments, 
and propose an experimentally feasible scheme to directly 
detect the observable effect induced by the noncommu- 
tative property of the non-Abelian gauge potentials in 
a trapped-ion system^. We show that a true non- 
Abelian gauge potential naturally arises in connection 
with the geometric phase associated with two degener- 
ated dark states in a four-state atomic system interact- 
ing with three pulsed laser fields. The system under con- 
sideration is a four-level atom interacting with resonant 
external laser pulses. Under an appropriate configura- 
tion, there exist two-fold degenerated dark states and 
then a non-Abelian gauge potential naturally arises. We 
design two specific closed loops C\ and C2 in the param- 
eter space, which generate two different geometric phase 
factors Ui and U2, respectively. Consider the composite 
path in which one traverses first C1 and then C%. After 
a cycle, a phase factor of U = U2U1 is generated. In 
contrast, the other phase factor of U' = U1II2 is accu- 
mulated if one traverses first Ci and then C\. We find 
that the population of the atoms in a specific state at the 
end of the two evolutions U and U' can be significantly 
different. For instance, the population difference can be 
as high as 43% for the typical experimental parameters. 
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Note that the observe precision of an atomic state can be 
higher than 99%, and thus the proposed noncommutative 
effect can be detected unambiguously. Since the varia- 
tions of the parameters associated with the non- Abelian 
geometric phases can also be easily manipulated via the 
laser beams, the proposed scheme is very promising for 
detecting the fundamental non- Abelian gauge structures. 

We start by reviewing a general framework of non- 
Abelian adiabatic geometric phase 0, [|[ and its possible 
detection in nuclear quadrupole resonance 16]. Consider 
a family of Hamiltonians H(x^) (m = 1,2, ■■■ , n) de- 
pending continuously on parameters Xya &U of which have 
a set of N degenerate levels. Let r] a be an A^-fold degen- 
erate set of orthonormal instantaneous eigenstates of the 
Hamiltonian H described by H \q a ) = E\rj a ) with E being 
the eigenvalue. Under the adiabatic condition the instan- 
taneous state \tp a ) will not overflow the state vector space 
spanned by rj a - So \ip a ) can always be expanded as a su- 
perposition of |?y a ): \ip a ) = T,b\Vb)Uba- Substituting the 
wave function \ip a ) into the time-dependent Schrodinger 
equation ij^\ij) a ) = H\ip a ), we have 



U ba 



(r/b\i)c)U c 



(1) 



where a trivial overall dynamical phase term has been 
dropped. A gauge structure appears when the Hamilto- 
nian H varies as the parameters x^(t) varv slowly with 
time t. We can define a gauge potential 



d 

and have (i] a \i] b ) = Y, fl A abfl (dx' 1 /dt). 
(2), it is straightforward to obtain 



(2) 



Integrating Eq. 



Uab = 



Vexp - / A^dx 



(3) 



ab 



where V denotes the path-ordered operator. The quan- 
tity A denned in Eq.© indeed plays the role of a gauge 
potential and the gauge group corresponds to the unitary 
freedom in choosing the basis states r\ a [2J ■ An intriguing 
feature of the gauge potential A lies in that it depends 
only on the geometry of executed path in the space of 
degenerate levels. 

In an earlier experiment, Tycko reported the effect 
of geometric phase on the nuclear-quadrupole-resonancc 
spectra [lf|. The effective Hamiltonian describing Ty- 
cko's experiment is given by 

H = (S-B) 2 

= (5a, sin 9 cos tp + S y sin 9 sin <p + S z cos 9) 2 B 2 ,(4) 

where S denotes the nuclear spin operator for the sample 
used in the experiment (a spin-| CI atom in a NaClO^ 
crystal is considered) , and B represents an external mag- 
netic field. There always exist a pairwise degeneracy of 
states in the parameter space because the Hamiltonian is 



invariant under the operation S — > — S, thus non- Abelian 
gauge structure may appear in the system. To address 
clearly the gauge structure, one can parameterize the 
Hamiltonian (0| as follows 



H = B A e 



-i9Sy g2 giOSy ^i<pS z 



(5) 



As a result, the instantaneous eigenstates can be directly 
written as 



Wa. 



_ e -iipSz e -i6S v 



(0) 



where |o) is one of the atomic state in the set {| ± |), | ± 
i)}. Here the state \m) represents the eigenstates of the 
spin operator S z defined by S z \m) = m\m) and the states 
| ± to) form a doubly degenerate sector. Substituting 
Eq.(J5J) into Eq.©, the corresponding gauge potentials 
are given by 



Aabip = (-i)(a\(cos9S z 
Aabe = (-i)(a\S y \b). 



sin9S x )\b), 



(7) 



As clearly shown in Ref.[3j, an Abelian structure occurs 
for \m\ = | sector although those states are doubly de- 
generate, while a true non- Abelian structure may appear 
for the sector \m\ — \. Explicitly, the gauge potentials 
for the later case is given by, 



A = —i 



1 „ a 

-a z cost) - —o x sm( 



dip 



-a r d9 



, (8) 



where a = S + 1/2 and cr XiV ^ z are Pauli matrixes. 

To disclose unambiguously the non- Abelian character- 
istic of the geometric phase, we must detect the physi- 
cally observable effects induced by the noncommutativity 
of the gauge structure. To this end, here we consider two 
closed loops C\ and C%, both starting and ending at the 
same point xo m the parameter space, where the two 
corresponding evolution operators are denoted as U± and 
U2, respectively. For the composite path that the sys- 
tem first traverses along loop C\ and then C2, the total 
phase factor of U — U2U1 is generated. If the counter- 
ordered evolution is executed, i.e., an evolution with first 
loop Gi and then C\ , the phase factor generated is then 
given by that of U = XJ-JJ-^- 111 a general experiment, if 
the observable effects are different for the phase factors 
of U and U , then Non- Abelian geometric structures are 
confirmed there. Otherwise, only the Abelian part of the 
gauge structure is observed even if the quantum states of 
the system are degenerate. 

From the above discussion, it is easy to understand 
that only the Abelian part of the gauge structures has 
been experimentally observed in the Tycko's experiment, 
as has been addressed in Ref. [3]. In the Tycko's ex- 
periment, the angle 9 is fixed as a constant (cos6> = 
l/v3). It is due to that only the rotation around a 
fixed axis can be implemented fast enough to gener- 
ate observable phase shift in NMR systems. In this 
case gauge structure A is proportional to a fixed ma- 
trix cos 9^- — a sin 9^- and the evolution operator can be 
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written as U — exp[— i(p(a z /\/3 — cr x )]. Thus the U and 
gauge potential A in different time sequences are always 
commutable. From this point of view, we say that the 
gauge potential loses its non-Abelian characteristic and 
is of Abelianization. In order to demonstrate the non- 
Abelian characteristic in nuclear quadrupole resonance, 
the system should traverse a non-Abelian path in which 
both 9 and tp vary with time. Unfortunately, it is difficult 
in an NMR setup. 




FIG. 1: Schematic map of the four-level system interacting 
with resonant laser fields. States |0) and \i) (i = 1,2,3) are 
coupled by a resonant external laser pulse with effective Rabi 
frequency fli. 



Recently, with the rapid experimental progress in 
trapped ions [H Qj, ES 0, EH, lH HI , it is quite pos- 
sible to experimentally observe the non-Abelian charac- 
teristic of gauge potentials with current technology. We 
now start to present a scheme to detect the non-Abelian 
characteristic of the gauge potentials in such system. The 
system under consideration is an ion (or a set of ions) con- 
fined in a (linear) Paul trap. Each ion has three ground 
(or metastable) states |1), |2), |3), and one excited state 
|0). The ground states could be different hyperfine lev- 
els or in the same manifold but with different Zeeman 
sublevels, and they are coupled to the excited state |0) 
by resonant classical lasers with Rabi frequencies f2i, fl 2 
and f^3, respectively, as shown in Fig.l[12|]. The effective 
Hamiltonian for each ion in the rotating frame reads 

Hj = -?i(O 1 |0)(l| + fi 2 |0)(2| + n 3 |0)<3|) + H.c. (9) 

In a general case the Rabi frequencies fi^ can be pa- 
rameterized with angle and phase variables according 



to fii = fl sin 9 cos (j>e llpi , tt 2 
Q 3 = flcos8e iip3 , where fl = 



tt sin 6 sin Se 1 ^ , and 



Two adiabatic eigenstates of Hi with zero eigenvalue 
(dark state) are then found to be 



sin</>e 4V |l) -cos^e iv |2), 



\D 2 ) = cos6»cos 9 !.e' i¥ '|l)+cos6lsin(/.e i¥ '|2)-sin6'|3), (10) 

where we have assumed that the laser phases satisfy the 
relations <p = tp 3 — tp% = tp 3 — tp%. Substituting Eq.(10) 



into the formula (2), we have, 
A = 0, 

o - cos e 



a 4 



COS ( 







i 
i cos 2 9 



(11) 



Then it is straightforward to find that the gauge poten- 
tials are given by 



A = S„A, l d X p 



-I 



a z I dip — ia y cos Odcj^Vl) 



where / is a 2 x 2 unit matrix. Similar to the above 
analysis in NMR system, the non-Abelian characteristic 
occurs only when both the angle <f> and phase ip vary with 
time. For instance, if the phase tp is fixed, we have 



U = exp(z7tr y ) 



cos 7 sin 7 
— sin 7 cos 7 



where 7 = J cos 9d<f>. This evolution and its geometric 
phase have been adressed in Ref. (2(| . However, it is clear 
that the non-Abelian characteristic (noncommutativity) 
of the potential A is lost in this specific case. There- 
fore, rather than non-Abelian geometric phase factor, the 
scheme proposed in [26] can only be used to observe the 
Abelian part of the geometric phase, although the two- 
fold dark states remain degenerate as the Hamiltonian 
varies with time. 

To observe the true non-Abelian structure, we propose 
two specific non-Abelian closed loops G\ and C2 in which 
both the phase tp and angle <fi vary with time. In the 
closed loop Ci , we assume that the Rabi frequencies are 
given by 

fii = fio/(*)> 

n 2 = n f 2 (t), 
n 3 = o e"* 2/T V v , 

where ip = ^t, and f(t) is set as 



(13) 



f{t) = I COBfg), 



-T<t<T, 

otherwise. 



(14) 



In this case, the parameters 9{t) and <j)(t) are determined 
by 

in 2 | firt\ 
tan^) = ^-=cosf-j, 



tan 9{t) 



|fil| 2 + |fi 2 | 2 _ COS 2 (f)+COs4(f) 



exp(-^) 



(15) 



From time — r — * t, the parameters (9, <p, p) vary from 
(0, 0,— 7r) to (0,0, 7r), and it just accomplishes a closed 
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path in the parameter space. It is remarkable that all 
parameters (9, 4>, ip) are completely independent of the 
laser parameters and r, provided that the adiabatic 
approximation is reasonable, thus it is easy to design the 
required laser pulses. By substituting Eq. (fl5|) into Eq. 
(fl~2|) . the gauge potential Ai generated by the loop C\ 
can be derived, and the associated evolution operator is 
given by 



U x = Vexp \- j( AZdxfj 



(16) 



We set the Rabi frequencies for the loop C2 as follows, 



fii = no/(t), 

3 = Qoe^-^^e^, 



(17) 



where a and j3 are two variables to differentiate the loops 
C\ and C2, and (3 can also be considered as a time delay 
among pulses ^3 and Ox, ^2- In this loop both the start 
and end points in the parameter space (#, 0, ip) are the 
same as those in the loop C\. Similarly, a gauge potential 
A2 will be generated in this loop C2, and the associated 
evolution operator is given by 



U 2 = Vexp 



J Co 



(18) 



The analytical result of integral in Eq. (fTB"]) (or Eq.( 
[15]) ) is difficult in a general case, because each segment 
in the time-ordered exponential loop integral does not 
commute with the next. Nevertheless, the integral can 
be evaluated easily by numerical methods. 

We now start to show that the composite path formed 
by the loops C\ and C2 has noncommutative feature. 
It is straightforward to check that the initial eigenstates 
|^i,2)i and the final eigenstates \Di^) f for both pathes 
C\ and C2 are given by \D\)i = l-Di)/ = — 12) and 
\D2)i — I-D2)/ = We may initially prepare the sys- 
tem as = \D2)i- Then we let the effective Hamil- 
tonian undergo a composite path, first C\ and then C2. 
The total evolution operator is a two-by-two matrix de- 
noted as U = U2U1. After this evolution, the final state 
becomes 

I*) / - C/12IA) / + U 22 \D 2 ) f = -U 12 \2) + U 22 \l), (19) 

where U\ 2: U 22 are the elements of the matrix U. On the 
other hand, if we let the Hamiltonian tracks the counter- 
ordered path, i.e., first C 2 and then Ci, the final state 
becomes 



I* >/ = ~U 12 \2) 



033 1 1) 



(20) 



where U 12 and U 22 are the elements of U = UiU 2 . De- 
note P = IC/22I 2 and P = \U 22 \ 2 , which correspond to 
the possibilities to find the final state in the state |1) after 



composite paths C 2 C± and C\C 2 , respectively, the pop- 
ulation difference Pd between the two composite paths is 
given by 



Pi. = P' - P = \U 22 \ 2 - \U: 



>22\ 



(21) 



The most essential feature of non-Abelian gauge struc- 
ture (the noncommutativity character) can be unambigu- 
ously observed if the population difference Pd 7^ 0. 

To show clearly the observable effect induced by the 
noncommutativity feature of the non-Abelian geometric 
phases, we here numerically calculate the population dif- 
ference described by Eq. pTj) for some typical parameters. 
The populations P, P', and the difference Pd between 
them as functions of the parameters a and /3 are plot- 
ted in Fig. 2 (a) and (b). In the calculations for the 
evolution U (U 1 ), we let U± (U2) implement during time 
— r — > r, while U2 (Ui) implement from r to 3r. The ef- 
fects induced by the noncommutative property are clear 
since the difference Pd can be non-zero. To have a quan- 
titative idea about the difference between the evolution 
U and U', the population difference Pd as a function of 
P for a = 7 is shown in Fig. 2(c). In this case Pd — 43.2% 
is derived for = 0.5. 

Nowadays the measurement technique of population 
of atomic states has been well developed and the de- 
tection precision can be higher than 99%. This high- 
precision measurement really makes it experimentally 
feasible to directly "see" the noncommutativity of non- 
Abelian gauge potentials by measuring the populations 
of atomic states. 




FIG. 2: (Color online) (a) The populations of atomic state 
|1) versus variables a and /3. P (P') is the population of the 
state |1) at the end of the evolution U (U ). (b) The difference 
P' — P of the two populations shown in (a) as a functions of 
the parameters a and 0. (c) The populations P, P' and P' —p 
versus the parameter (3 when a — 7. 

We now turn to demonstrate that the adiabatic condi- 
tion crucially required in the proposed scheme is well sat- 
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FIG. 3: (Color online) (a) The populations Pdi (Pd2) for the 
dark state \D\) (I-D2)) during a complete composite evolu- 
tion U. (b) The populations Pbi {Pbi) for the bright state 
\Bi) (\B 2 )). (c) The maxima Pg ax of the sum P B i and P B2 
during the whole evolution time versus parameters /3. The 
parameters flo = 100 and r = 2. 



isfied. Besides the two degenerate dark states described 
in Eq. pi?)) . there exist two other bright states with eigen- 
values ±Q in the system described by the Hamiltonian 
([9|), and the corresponding eigenstates are given by , 

\B 1 )= -J^ [sin 61 cos ^ 1 1) + sin 8 sin cj)e iip \2) 
V2 

+ cos0|3) + e iv |O)], 

\B 2 ) = -Usmflcos^ll) + sin 6» sin ^12) 
V 2 

+ cos0|3)-e i *'|O)]. (22) 
The adiabatic approximation is well satisfied if the over- 



flow from the initial dark states to the two bright states 
(|22| can be negligible during the whole evolution. In view 
of this point, we numerically calculate the populations of 
all dark and bright states. The results for the parameters 
a = 7 and f3 — 0.5 are plotted in Fig. 3(a) and (b). It 
is shown that the populations for both bright states can 
be negligible. In particular, we calculate the maxima 
(P B vIax ) of the sum Pbi and Pb2 during the whole evolu- 
tion time versus the parameter (3 with the other param- 
eters being the same as those in Fig. 2(c). We note that 
the sum of the populations for the two bright states are 
always below 1.0% during the whole evolution. Therefore 
it clearly shown that the adiabatic approximation used in 
the calculations of data in Fig. 2 (c) is reasonable. Usu- 
ally the adiabatic approximation is well-satisfied when 

fv» i [liSzI- 

In summary, we have proposed an experimentally feasi- 
ble scheme to detect a true non-Abelian geometric phase 
effect in trapped ions coupled with laser beams. By de- 
signing two specific composite cyclic evolutions formed 
by U2U1 and U1U2, we show in detail how the effect in- 
duced by the noncommutativity of non-Abelian gauge 
structures can be observed through detecting the popu- 
lation of the internal states of trapped ions. The high- 
precision detection of the atomic states recently devel- 
oped in the field of quantum information opens the great 
possibility to directly "see" the noncommutative feature 
of the gauge potentials. It is quite promising for exper- 
imentalists to implement the interest idea presented in 
this work. 
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